Suppression of rectification at metal-Mott-insulator interfaces 
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Charge transport through metal-Mott-insulator interfaces is studied and compared with that 
through metal-band-insulator interfaces. For band insulators, rectification has been known to occur 
owing to a Schottky barrier, which is produced by the work-function difference. For Mott insulators, 
however, qualitatively different current-voltage characteristics are obtained. Theoretically, we use 
the one-dimensional Hubbard model for a Mott insulator and attach to it the tight-binding model 
for metallic electrodes. A Schottky barrier is introduced by a solution to the Poisson equation with a 
simplified density-potential relation. The current density is calculated by solving the time-dependent 
Schrodinger equation. We mainly use the time-dependent Hartree-Fock approximation, and also use 
exact many-electron wave functions on small systems for comparison. Rectification is found to be 
strongly suppressed even for large work- function differences. We show its close relationship with the 
fact that field-effect injections into one-dimensional Mott insulators are ambipolar. Experimentally, 
we fabricated asymmetric contacts on top of single crystals of quasi-one-dimensional organic Mott 
and band insulators. Rectification is strongly suppressed at an interface between metallic magnesium 
and Mott-insulating (BEDT-TTF)(F 2 TCNQ) [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene, 
F2TCNQ=2,5-difluorotetracyanoquinodimethane]. 



I. INTRODUCTION 

When a new function of some material like an organic 
molecular crystal is sought, one often considers a possible 
electronic device made of the material. Around the de- 
vice, there always exists an interface between two materi- 
als with different transport properties and different work 
functions. Charge transport through metal- insulator or 
metal-semiconductor interfaces has been studied from 
long ago* 1 * 2 - However, most of the theories are restricted 
to metal-band-insulator interfaces. 3 Now it has become 
clarified that many organic conductors and insulators 
are strongly correlated electron systems. 4 Their insulat- 
ing phases cannot be explained on the basis of band 
structures. Electron-electron interactions play an essen- 
tial role in the insulating phases of the (BEDT-TTF) 2 X 
salts. Then, conventional theories cannot be applied 
to charge transport through metal-Mott-insulator inter- 
faces. 

As is well known, two different materials generally have 
different work functions. When they are attached to each 
other, the energy levels are modified around the inter- 
face by transferring electrons to match their Fermi levels. 
This is often called band bending, and forms a Schottky 
barrier at the metal-insulator interface. The electronic 
state in the close vicinity of the interface can be modified 
physically and/or chemically and depends sensitively on 
the two materials^ Nevertheless, the overall band struc- 
ture is governed by the long-range Coulomb interaction, 
i.e., by the Poisson equation if the band structure can be 
continuously treated! 3 - 

When a Schottky barrier dominates the characteris- 
tics of a field-effect transistor fabricated on an organic 
insulator, the field-effect characteristics indeed depend 



on whether it is a band insulator or a Mott insula- 
tor. For instance, when such transistors are fabricated 
on semiconducting carbon nanotubes, which are band 
insulators, the field-effect characteristics are generally 
unipolar unless the work function of the metallic elec- 
trodes is matched with that of the carbon nanotube^i^ 
Now, fine control of p-, n-, and ambipolar- type oper- 
ations is available in single-crystal organic field-effect 
transistors with the use of chemically tunable Fermi en- 
ergies in tetrathiafulvalene-tetracyanoquinodimethane- 
based organic metal electrodes*^ On the other hand, 
when metal-insulator-semiconductor field-effect transis- 
tor device structures are based on organic crystals 
of the quasi-one-dimensional Mott insulator (BEDT- 
TTF)(F 2 TCNQ), the field-effect characteristics are am- 
bipolar even if its work function is quite different from 
that of the metallic electrodes^ We have theoretically 
shown that backward scatterings at interfaces and umk- 
lapp scatterings inside the Mott insulator are balanced, 
leading to collective charge transport insensitive to the 
band bending. 12 

Then, how about a simpler structure of just a metal- 
insulator interface without gate electrode? In order for 
the current-voltage characteristics to be observed, two 
different metallic electrodes need to be attached to an 
insulator in such a way that only one interface has a 
large work-function difference. The magnitude of the 
current density depends on the sign of the applied volt- 
age: the forward (reverse) voltage lowers (raises) the 
Schottky barrier at this interface, leading to a larger 
(smaller) current density for band insulators including 
ordinary semiconductors; 3 If this asymmetry is large, it 
can be regarded as rectification. Such rectifying nature 
is observed in single-crystal organic devices without gate 
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FIG. 1: (Color online) One-dimensional model for an insula- 
tor, to which two metallic electrodes are attached. 



dielectric layers.— Then, what happens to the current- 
voltage characteristics if they are replaced by Mott in- 
sulators? In this paper, we show that the rectification 
is strongly suppressed for Mott insulators. Preliminary 
theoretical results are presented in Ref. [l3|. 



II. ONE-DIMENSIONAL MODELS FOR 
METAL-INSULATOR INTERFACES 

We use the one-dimensional Hubbard model for a Mott 
insulator (the one-dimensional tight-binding model with 
alternating transfer integrals for a band insulator) at- 
tached to two different metallic electrodes represented 
by the one-dimensional tight-binding models with a com- 
mon transfer integral for simplicity. The total number of 
electrons is set the same as the number of sites. 

i i 
- X [*i,i+l(*) c L c i+l,<r +*i+l,t(*) c i+l,<7 < V >(!) 

where cj CT (cj jCr ) creates (annihilates) an electron with 

spin a at site i, n ia = c] CT Cj i<r , and rij = Yla n iv The 
site energy ej is set at 0l in the left electrode at 1 < 
i < (L e - l)/2, at in the insulator at (L e + l)/2 < 
i < L — (L c + l)/2, and at </>r in the right electrode at 
L — (L c — l)/2 < i < L, where L is the total number of 
sites and L c the number of sites in the electrodes (Fig.[T]). 
The absolute value of the transfer integral | ti^ + i{t) | is 
set at t c if either i or i + 1 is in the insulator and at 
t e otherwise. When we consider a Mott insulator, the 
on-site repulsion Ui is set at U in the insulator and at 
in the electrodes. For a band insulator, U =0 and the 
transfer integral t c above is replaced by t c — (—lySt. 

A scalar potential Vi is introduced above to account 
for redistribution of electrons at interfaces to form barri- 
ers, compensating work-function differences </>l and 0r in 
equilibrium. The periodic boundary condition is imposed 
on Vi. vq = vl- The applied voltage V is so denned that 
it is positive when the right electrode has a lower poten- 
tial (for electrons) than the left and the current (without 
multiplication of charge) flows to the right. In order to 
maintain the periodic boundary condition for finite V, 
we introduce the Peierls phase into the transfer integral, 



U,t+i{t) = t* +1 At) =| t ht+1 (t) | exp 



_i£5 J dt'E(t') 



where t denotes time, e the absolute value of the elec- 
tronic charge, a the lattice constant, and E(t) the aver- 
aged electric field denned by E(t) = — V/(La)M 

By adding the vector potential introduced into the 
Peierls phase to the scalar potential Uj, the total poten- 
tial ipi is given by ipi = Vi — V(i/L — 1/2). Although 
the potential ipi is defined on lattice points, we solve the 
Poisson equation in the continuum space, 



dx 2 



■Vp ((n) - I) for L c /2 < x < L - L c /2 , (3) 



and d 2 ^/dx 2 =0 otherwise, where the potential ip and 
the expectation value of the electron density per site (n) 
are functions of x, and Vp comes from the long-range 
Coulomb interaction. In order to match the Fermi levels, 
we set the boundary condition, i.e., the potentials in the 
metallic electrodes, as 

i>{x) = -<f) L + V/2 for < x < L c /2 , 

i/j(x) = -0r - V/2 for L - L /2 < x < L . (4) 

In order to solve the Poisson equation analytically, we 
introduce a simplified density-potential relation, which 
is known to work well for the potential distribution even 
in strongly correlated electron systems^ 



dn(ip) 
dip 



W 



(5) 



with a constant compressibility n, the bandwidth W, 
and the gap A. We have confirmed, by using the ex- 
pectation value (n) with respect to the wave function 
obtained from the time-dependent Schrodinger equation 
and by solving the Poisson equation for ip simultaneously, 
that the current-voltage characteristics are qualitatively 
unchanged (Compare the present results with those in 
Ref. [H). For the expression of ip(x), see Appendix. 

First, we set V at zero and obtain the ground state ei- 
ther exactly (for small systems with L <14) or in the un- 
restricted Hartree-Fock approximation. Then, we set V 
at a finite value to solve the time-dependent Schrodinger 
equation. For exact many-electron wave functions, the 
exponential evolution operator with time slice dt=l0~ 2 is 
expanded to the 15th order. For self-consistent Hartree- 
Fock wave functions, the evolution operator with time 
slice dt—10~ 4 is decomposed with the help of the Suzuki- 
Trotter formula so as to be accurate to the order of dt 2 . 
The current density J(t) is averaged over the period, 
< t < At with At = 2TrhL/(4eV)^- to give I, where 
J(t) is given byi^^ 



A*) = (1/£)X [^+M(*)4+i,a<V -**M+i(*)4« 



III. NUMERICAL RESULTS 



Cj+1,0 



(6) 



(2) 



The analytic solution ipi given in Appendix plus the 
site energy ej — Ui/2 is plotted in Fig. O The energy 
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FIG. 2: (Color online) Energy levels, ti + ipi ± t/j/2, where 
ipi is the solution to the Poisson equation © with the simpli- 
fied density-potential relation ((5J) for L=100, L =49, W=4.0, 
A=0.33 (corresponding to (7=1.6, <5t=0 or (7=0, <ft=0.0825 
with t c =l), Vp=0.05, L =O.65, and <A R = -0.165. A Schottky 
barrier is formed at the left interface, and a positive (negative) 
V for right- (left-)going electrons corresponds to a forward 
(reverse) voltage lowering (increasing) the barrier height. 



levels for doubly occupied sites, + ipi + Ui/2 in the case 
of a Mott insulator, are also plotted as a guide. This 
result is used later in calculating the time evolution of the 
system. Here, the work function of the right electrode is 
set to be matched with the top of the lower Hubbard or 
valence band. The work function of the left electrode is 
set to be much lower than that of the insulator or that of 
the right electrode. A Schottky barrier is formed at the 
left interface. Its barrier height is lowered (increased) for 
right- (left-)going electrons, so that a positive (negative) 

V is a forward (reverse) voltage. 

In this situation, one usually expects that the absolute 
value of the current density /, | I |, is large (small) for 

V >0 (V <0), i.e., rectification behavior.— It is indeed 
reproduced for a band insulator, as shown in Fig. [3fb). 
However, the current-voltage characteristics for a Mott 
insulator are qualitatively different, as shown in Fig.[3ja). 
Namely, the absolute value of the current density is in- 
sensitive to the sign of the voltage, so that the current 
density I is almost an odd function of V for a wide range 
of (pi,. The rectifying action is strongly suppressed for 
the Mott insulator. This property is maintained even for 
large 4>l- We have performed numerical calculations with 
different sizes L, L c and with different Coulomb param- 
eters Vp to find that this property is quite robust. 

In order to confirm that this property is not an ar- 
tifact of the approximation, we then use exact many- 
electron wave functions on small systems and compare 
their current-voltage characteristics with those by the 
time-dependent Hartree-Fock approximation in Fig. [4] 
Here a large Coulomb parameter Vp is used to accommo- 
date the band bending within the small insulator of size 
L — L c . We also use U—Q and St=0.375 for a band insula- 
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FIG. 3: (Color online) Current-voltage characteristics for (a) 
Mott insulator with (7=1.6, 8t=Q, and (b) band insulator with 
(7=0, 5t=0.0825, obtained by the time-dependent Hartree- 
Fock approximation for L=100 and L c =49. Other parameters 
are t e =t c =l, Vp=0.05, <Ar,= —0.165, with different <j>L values 
as indicated. Rectification that is clearly seen for the band 
insulator is suppressed for the Mott insulator. 



tor with gap A=1.5, which is close to the gap in Fig.[4jb), 
to confirm that rectification survives in the band insula- 
tor on this small system. In contrast, rectification turns 
out to be suppressed in the Mott insulator irrespective of 
whether exact many-electron wave functions are used or 
the time-dependent Hartree-Fock approximation is em- 
ployed, although they show a quantitative difference. 

The similar qualitative difference between Mott and 
band insulators is demonstrated by field-effect carrier 



injections 



11.12 



The field-effect characteristics are always 



ambipolar for Mott insulators and unipolar for band in- 
sulators when the work functions are different between 
the insulator and the metallic electrode. This field- 
effect property (with drain current Ip> , drain voltage Vp, , 
and gate voltage Uq in the notations of Ref. Qjl) can 
be shown to be closely related to the present interfa- 
cial property. The present result of / as a function 
of V is denoted by I = f(V). Then, in the geome- 
try of Ref. [l2|, the left (right) interface tends to pro- 
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FIG. 4: (Color online) Current-voltage characteristics for 
(a) exact many-electron wave function, and (b) Hartree-Fock 
wave function, on a small system with L=14 and L c =5. Other 
parameters are (7=2.7, 8t=0, t e =t c =l, Vp=5, 0l=O.75, and 
4>-r~ —0.75. The overall characteristics showing suppressed 
rectification are obtained for both wave functions. 



duce a forward (backward) current density I\ (J 2 ) by the 
potential difference U G + V D /2 (U G - Vjj/2). The to- 
tal current is approximately given by In oc I± — J 2 = 
f(U G + Vd/2) - f(U G ~ Vd/2) ~ V D /' (C/g) for small Vb, 
where f'(V) is the derivative of /(V). Therefore, the sup- 
pressed rectification [/(V): odd function of V] leads to 
ambipolar field-effect characteristics \f'{U G ): even func- 
tion of U G \. The above relation between the interfacial 
property and the field-effect property through differenti- 
ation approximately holds for band insulators also. 

For the field-effect characteristics, the gate-bias polar- 
ity possessing the higher Schottky barrier at the metal- 
Mott-insulator interface is accompanied by more devia- 
tion from half filling (i.e., weakened umklapp scattering) 
inside the channel than that for the opposite gate-bias 
polarity^ Such a counterbalance is observed in detail 
and found in a very wide parameter space spanned by 
the work-function difference, the bandwidth difference, 
and the Coulomb parameter in the Poisson equation. In 
the present case, a very similar mechanism works: the 



bias-polarity dependence of the present charge-density 
distribution around the Schottky barrier is quite simi- 
lar to the corresponding one of the field-effect transistor. 
Charge transport through a Mott insulator is not simply 
governed by the interface but determined as a whole in- 
cluding the interfacial and bulk regions. It is therefore 
much less sensitive to the details of the interfacial barrier 
potential than charge transport through a band insula- 
tor. Whether the voltage is forward or reverse is insignif- 
icant at metal-Mott-insulator interfaces. Such collective 
charge transport turns out to be realized by balancing 
the barrier effect with the correlation effect, as shown in 
Fig. 9 of Ref. GJ. 



IV. EXPERIMENTAL RESULTS 

In the experiment, we prepared single crystals of 
two kinds of quasi-one-dimensional organic charge- 
transfer complexes; (BEDT-TTF)(F 2 TCNQ) and K- 
TCNQ (TCNQ=tetracyanoquinodimethane). The latter 
is known as a Peierls-type (band) insulator composed of 
segregated stacks of TCNQ anion radicals. The stack of 
the compound is strongly dimerized along the stacking 
axis at room temperature^ The former is a Mott in- 
sulator composed of side-by- side arrangement of BEDT- 
TTF cation radicals. The molecular arrangement is free 
of dimerization down to low temperature (~ 4 K). We 
fabricated asymmetric contacts on top of the single crys- 
tals by vacuum deposition with a gap of about 100 /im, 
to measure the current along the quasi-one-dimensional 
stacks or chains. Cathodes are fabricated with 100 nm of 
magnesium (work-function </>= 3.66 eV), which is coated 
with thin (3 nm) silver layers to prevent the oxidation of 
magnesium. The contacts form Schottky junctions with 
both crystals. In contrast, anodes are fabricated with 
50 nm of gold (<f>= 5.1 eV) for (BEDT-TTF) (F 2 TCNQ) 
crystals and 50 nm of silver (0= 4.26 eV) for K-TCNQ 
crystals. The contacts form ohmic junctions with the re- 
spective crystals at room temperature. The DC current- 
voltage characteristics were measured with the use of the 
semiconductor parameter analyzer (Agilent E5270). 

Current-voltage characteristics of the devices are 
shown in Fig. [5] It is found that the device with K- 
TCNQ exhibits rectifying nature with the rectification 
ratio of about 5 [Fig. Ob)]. The results clearly show the 
formation of typical Schottky junctions in K-TCNQ with 
magnesium. In sharp contrast, the device with (BEDT- 
TTF) (F 2 TCNQ) does not show rectification, while it ex- 
hibits distinct nonlinear features in the low- voltage range 
at both polarities [Fig. 02a)]. The observed nonlinearity 
should be ascribed to the cathode characteristics with 
magnesium, since the device with gold contacts for both 
the cathode and the anode show much more conduc- 
tive characteristics. It should be noted here that elec- 
tron correlations are actually strong in K-TCNQ, but the 
dimerization-induced backward scattering does not allow 
the balance between the barrier effect and the correlation 
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used for a Mott insulator and the tight-binding model 
with alternating transfer integrals for a band insula- 
tor. For small systems, we have compared numerical 
results obtained through the time-dependent Hartree- 
Fock approximation with those of exact many-electron 
wave functions. The qualitative characteristics regard- 
ing the current-voltage asymmetry are not lost by this 
approximation, which is used for large systems. Rec- 
tification is shown to be strongly suppressed at metal- 
Mott-insulator interfaces even with large work-function 
differences. This property dominated by an interface is 
shown to be closely related to the previously explained, 
ambipolar field-effect characteristics of a one-dimensional 
Mott insulator with two interfaces. This collective charge 
transport through metal-Mott-insulator interfaces is a 
consequence of the fact that backward scatterings at in- 
terfaces and umklapp scatterings inside the Mott insula- 
tor are balanced. The suppression of rectification is ex- 
perimentally observed in the current- voltage characteris- 
tics of the quasi-one-dimensional Mott insulator (BEDT- 
TTF)(F 2 TCNQ), with magnesium cathode forming a 
Schottky junction and gold anode forming an ohmic junc- 
tion, at room temperature. 
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APPENDIX: POTENTIAL ip{x) 



FIG. 5: (Color online) (a) Current-voltage characteristics of 
(BEDT-TTF)(F 2 TCNQ) with magnesium cathode and gold 
anode at 295K. (b) Current-voltage characteristics of K- 
TCNQ with magnesium cathode and silver anode at 250K. 



effect, leading to the band-insulator-like behavior. Thus, 
we found that these experimental results are consistent 
with the theoretical arguments presented above. 



V. SUMMARY 

In order to study the effect of electron correlation on 
charge transport through metal-insulator interfaces, we 
consider a one-dimensional Mott insulator, to which two 
metallic electrodes with different work functions are at- 
tached, and observe its current-voltage characteristics. 
The spatial dependence of the potential including Schot- 
tky barriers at interfaces is theoretically modeled by an 
analytic solution to the Poisson equation with a simpli- 
fied density-potential relation. The Hubbard model is 



We adopt the simplified density-potential relation 
which implies 

n(ip) = 2 for - oo < ip < -W/2 , 

n(ip) = 1 — n(tp + A/2) for - W/2 < ip < -A/2 , 

n(ip) = I for - A/2 < V < A/2 , 

n(V>) = 1 — ft (ip — A/2) for A/2 < V < W/2 , 

ra(V>) = for W/2 < tp < oo . (A.l) 

Then, the solution to the Poisson equation §5§ is given 
by 



ipix) = —A/2 ± Co sinh \J k Vp (x — xq) 
with sign + for x < xq and — for xq < x, and 



for I < n < 2 
(A.2) 



ip(x) — A/2 ± cq sinh \J nVp(x — xq) 



for < n < 1 



(A.3) 

with sign + for xq < x and — for x < xq, where Co and 
Xq are constants. Taking the boundary condition ((?]) into 
account, we obtain (with definitions 0l(F) = </>l — V/2 



6 



and <Pr(V) = </>r , + V/2) for example: 

i) if - y < -^l(VO < y and - ^ < -<f> R (V) < j 



where Xq is the solution to 
^kV p (x - L e /2) 



-0 L (F)+A/2 



ip{x) = 



MV){x -L + Le/2) - 4>k{V)(x - L e /2) 



L-L c 



for L c /2 < x < L - L c /2 



(A.4) 



sinh [^/kVp(L - L c /2 - x )] -A/2 + MV) ' 

(A.9) 

and Co is given by 



co 



W A A A 

ii) if - — < -MV) < -- and - - < -MV) < - 



-fa(y) + A/2 
/^Vp(x - Le/2) 



(A.10) 



ip(x) = —A/2 + Co sinh \J nVp (x — xq) 

for L c /2 < x < xq , 

v (A/2)(x - L + L c /2) - R (^)(x - xq) 

Vw = ; — rn 

L - Le/2 - Xq 

for So < x < L — L /2 , 

where Xo is the solution to 

sinh [V^Vp{x q - Le/2)] _ -A/2 + L (VQ 
% /^Pp(L - Le/2 - x ) ~ -MV) + A/2 

and Co is given by 

~MV) + A/2 



W A A W 

iv) if - — < -MV) < -j and - < -MV) < Y 



(A.5) 



tp(x) = —A/2 + Co sinh y «Vp(x — Xo) 
for L e /2 < x < xo , 



A 



for xq < x < xi , 



(A.6) 



ip{x) = A/2 + Co sinh \J nVp (x — x\) 

for xi < x < L — L e /2 , 
where xo and X! are the solutions to 



sinh [y/nVp{x - L c /2)] _ -A/2 + M V ) 



Co 



V^(L - L c /2 - x ) 



(A.7) 



VkVp(xi - x ) 



A 



(All) 
(A.12) 



and 



A AW A 

hi) if - - < -MV) < y and - — < ~MV) < -y , 



VkVp(xi - x ) 



V»(x) 



-( A/2) (x- Le/2) + ( /) L (^)(x-x ) 
x - L e /2 



sinh [V*(L - L c /2 - x x )] -<^r(F) - A/2 ' 

(A.13) 

and cq is given by 



for L c /2 < x < xo , 

V>(x) = —A/2 — Co sinh \J nVp(x — xo) 
for xq < x < L — L /2 , 



Co 



/ kVp(xi - x ) 



(A.14) 



(A. 
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